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parts behind of the enclosure simulated based on numerical methods, FDTD, 
MOM and equivalent circuit for better analysis of electromagnetic 
Keywords: interference. After introducing the proposed structure and presenting the 
curvature theory, simulation results are displayed and compared in the 


Apertures selected frequency range for three numerical methods. It has been shown that 
Complex structure increasing the number of apertures by reducing the size, increases the 
Electromagnetic interference effectiveness of the protective shield. However, increasing the number of 
Shielding effectiveness resonances by increasing the apertures indicates the importance of studying 


the equipment more precisely before choosing the structure of enclosure. 
We present a complex structure for the enclosure and the different number 
and dimensions of apertures with different materials were investigated for 
analyzing the effect of shielding on electromagnetic interference. The 
necessity of choosing a more effective enclosure according to the frequency 
of the equipment is specified. Finally, three methods of numerical solution, 
FDTD, MOM and circuit equal comparition were performed with measured 
value. Changes in the Shielding effectiveness and the number of resonant in 
the frequency range were determined. The exact examination of equipment 
requires shielding and their frequency and the type of inside-to-outside 
communication device before choosing shieldin is important. We used a 
comparison of three numerical solution methods for examining the field 
distribution in a complex structure enclosure with different apertures and 
different materials. In the majority of cases, the proximity of the measured 
values in this frequency range with the MOM curves shows the performance 
of this method in complex structures. 
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1, INTRODUCTION 

At present, the use of metal enclosures has increased for devices operating at high frequency [1]. 
In the frequency of operation, if the wavelength range and the dimensions of the enclosure are close, 
the effect of the enclosure in the interference should be well investigated. Electromagnetic interference (EI) is 
a complex, continuous and random signal, and requires accurate analysis and measurement [2]. In fact, 
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the presence of apertures in a metal enclosure is due to the presence of connectors and communication cables 
in the shielding enclosure. Different methods of numerical solution are used to compute Shielding 
effectiveness (SE) in complex compartments [3]. To calculate the SE in_ electromagnetism, 
numerical solution methods such as Method of Moments (MOM), Finite-difference time-domain (FDTD) 
method and Equivalent Circuit Methods (ECM) are used. Solving the integral equations in electromagnetic 
with using numerical methods is studied for more than 4 decades. Simple shapes in the form of square and 
rectangular cylinders, circles, etc. have been analyzed over decades [4]. The use of basic functions in many 
of these numerical methods and transforming the integral equation into a linear system and finally, solving it 
with direct methods and repetition is the ideal solution to achieve the desired goal [5]. The simplest method 
for solving electromagnetic problems is the FDM numerical solution method, which converts the differential 
equation into partial differential equations. In this way, the whole area is divided into nodes. The differential 
approximation of the original equation is written for them and solved with using algebraic algebraic 
equations matrix. This matrix is very large after the solution and does not support the open region, which is 
one of the important problems of this solution [6-7]. FDTD is a hybrid numerical method and used to solve 
electromagnetic problems. Because the rigorous numerical method and the asymptotic scheme do not work 
well, the combined method is more effective in solving problems in complex structures [8-10]. In the method 
of moments, as an advanced numerical solution, instead of using Maxwell differential equations, the Integral 
responses of Maxwell’s Equations are used. For this reason, If Green's function is known, the function of the 
MOM method is excellent for solving electromagnetic problems in complex structures [4] and [11]. 
To analyze the field distribution on the outer and inner surfaces of a metal enclosure, method of moment 
numerical method is an efficient method. In this method, by applying boundary conditions, the accuracy of 
field calculations has been improved [12]. An electromagnetic compatibility analysis of the components 
inside the enclosure is easy for simple structures with ordinary geometric shapes. However, the precise 
analysis of the performance of components in the complex structure with apertures is difficult and complex in 
terms of electromagnetic interference. So it needs to be reviewed and computed. With the advancement of 
technology, the growth of electronic equipment is faster than before [13]. Electromagnetic interference has 
become more and more complex in the environments, and as a result, the importance of electromagnetic 
compatibility analysis has increased. 

For this reason, the performance analysis of complex structures with different number of apertures 
and different shapes will provide a more accurate understanding of the destructive effects of electromagnetic 
interference. The enclosure with complex structure, with different dimensions of apertures, has been studied 
and compared with different numerical methods to investigate the interference effects of the equipment 
installed in the enclosure. The investigated shielding can create good conditions for accurate operation of 
high-frequency circuits in compressed and confined spaces [1, 2]. 


2. ANALYZING SHIELDING EFFECTIVENESS WITH DIFFERENT NUMERICAL METHODS 

At present, different shapes of metal enclosures are used to reduce electromagnetic interference. 
Shielding effectiveness is an important parameter for design and selection of the enclosure. This coefficient 
can be calculated by using equations | and 2 [12], [14]. 


SEg = 201log1) | dB (1) 
SEx = 20 logy a (2) 


2.1. Investigation of Electromagnetic Interference in Complex Curved Structure with Apertures by 
Equivalent Circuit Method 

The complex structure of the enclosure is shown in Figure |. Considering the appearance of these 
enclosures, it can be expected that conditions will be provided to use functional and compact enclosures 
instead of simple rectangular and square enclosures in the future. The internal curvature in both sides of the 
enclosures and the different number of apertures makes the analytical calculations complex and specific. 
The apertures act as the actual model for the arrival and departure of communication equipment in the 
enclosures. The effect of this structural form on the electric field and the intensity of the field are shown in 
Figure 2 [15]. The specification of the enclosure wall includes is, K the conductivity, 6 effect of the wall 
thickness, and €, the relative permeability coefficient. For exciting two short monopoles Antenna are used in 
the coordinates (x1, yi, 0) and (x2, y2, 0). L is length of the Antenna and r is the radius of the Antenna. 
The characteristics of the probes and Parameters for design of the enclosure are shown in Table 1. We use 
Coaxial cable for supply [16]. Metal shielding is one of the main tools for preventing electromagnetic 
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interference, which is a guarantee of electromagnetic compatibility in electronic equipment. SE coefficient 
specifies the shielding efficiency of a metal enclosure. The presence of apertures on the wall of enclosure 


changes the behavior of the shields. 





(a) (b) 


Figure 2. Distribution of magnetic (a) and electric (b) field in a complex special structure 


Table 1. Characteristics complex structure enclosure 


Caviti Current Probe 1 Current Probe 2 
L,=400 mm X,=90 mm X,=120 mm 
L,=500 mm yi=l110 mm y2=300 mm 
L,=300 mm L,=l mm L,=1 mm 
fo=100 mm 
E.= 22 

k = 10° s/m 


The simulation of the SE coefficient in complex structures with apertures is complicated by 
considering the solution of waveguide. Therefore, it is necessary to use the best method for calculating the 
SE coefficient in enclosure when resonance problems occur for large electric fields [3]. The evaluation and 
computation of field parameters at common boundaries between two different physical environments 1s 
necessary. The boundary conditions are shown in Figure 3. 





Figure 3. Boundary conditions in a complex structure 
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The tangential field components E are continuous along the whole boundary (Equation 3). The 
Vertical field components is calculated at the boundary of the two materials of Equation 9, where D is the 
density of the flux and p is the surface electrical charge density [18, 19]. 


V.D=p (3) 
VXE= (4) 
Eye = Exe (7) ) 
Din — Dan = Ps (=) (6) 
A =e = J, Eas, d; = Jq,D.4A, = Jz,eaV (7) 


The boundary conditions of B and H vectors analysis are necessary in environments with different 
physical properties. The static magnet is characterized by fundamental equations as Equations 8. 
Since divergence of B is zero in Equation 9, we conclude that the vertical component of B is continuous to 
crosses the boundaries. By the integral form of the Kearl equation, the boundary conditions of the tangential 
magnetic field component are calculated from Equation 10, which J,, is the surface flow density [18-20]. 


VxH=0 (3) 
Bin = Bon (9) 
Hie — Hoe = Jen (=) (10) 
hip=h"; =f Hds, dj=J,B.dA, j= J, Jad (11) 
Ch=—d+j, Ce=-<b, Sb=0, Sd=q (12) 


In the equivalent circuit method, the curved volume at the behind of the complex structure is excited 
by current probes. The equivalent impedance between these probes is obtained by using the power balance. 
Analysis of the Helmholtz equation for a curved surface with a volume V and a contact surface S has been 
performed. The electric Field is divided in two sections, a rotational part with a zero divergence (E,,) and a 
non-rotational Section with a zero Curl (F,,). Finally, the total magnetic field is computed and displayed. The 
total electric field equation is given (Equation 13). Since the non-rotational part is zero, we need to study the 
rotational section. This formula based on the Helmholtz equation is as follows (Equation 14). 
Given Equations (13) and (14), Equation (15) is obtained. Because E,and F,,are orthogonal, equations 16 
and 17 are calculated. Using the electric field expansion model (Equation 13) and using Equations 16 and 17, 
we can calculate the total impedance of the whole system (Zjj) as a two-port network with ports I and J, 
which is shown in Equation 19. 





E = Qo enEn + finFin (13) 
VxVXE-K°E = —jopj (14) 
Teg — ay — < J, En > (kn En — K? Eq) = —Jeopy “is 
n=0 (k2-k2)Ken ron nEn ie Ll 
Cn = Ga Kay SI En > (16) 
1 
li Te (17) 
— JOH yr00  SEnJj><En i> Lo <Fn J j><Fn Ji> 
Lij a Tl; dnb (k2-k2)KEn penis Din Kea (18) 
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2.2. Shielding effectiveness by FDTD Method 

Finite Difference Time Domain Method (FDTD) is an efficient method for solving Maxwell's 
equations. The main weakness of this method is the divergence error, which affects the accuracy of the 
numerical solution. The presence of stairs and curves in complex structures has a negative effect on the 
convergence of this method. To solve this problem, the combination of this method with the HO method 
(higher order FDTD) is used [15]. In Figure 4. the coordinates of the enclosure in the boundary region are 
visible. 









Zl (i,j + 1/2,k) 
G+1,j,kK) 


fi] G+1,j—1/2,k) 


Oe, 


Figure 4. Coordinates of the enclosure in the Boundary region 


It is assumed that the boundary specified with fi = (fi ,, iy, fi y) in Figure 4 is a normal unit vector. 
The components of the covariance of the electric and magnetic fields E™ and H™ on the boundary between 
the two regions A and B and for the values of u™ and e™ are calculated as follows. 


nx FSi KE nb. = 2 n.k (19) 
nXHGH=7n XH, e4fA.HA = e® fA. HB (20) 


If in the FDTD method, we assign the (i, j,k) = (iAu, jAv, kAw) character to the cell’s coordinates as 
in Figure 4. We assign a parameter to each cell with respect to its distance from the wall. This coefficient 
Br depends on the geometric shape of the boundaries. By applying the central finite difference and the 
above definitions, the discrete form of the Maxwell equations is as follows. These formulas are presented for 
H,, in two regions A and B. We extend H® by using Equation 21. All Maxwell equations are solved by 
numerical methods, and the analysis of electromagnetic interference by the SE is possible [15]. 


HAG K) = (1+ BA )Hu (i- 5k) + BAe Hu (i-5,i,k) (21) 
4HA(i,j,k) = Hy (ii-5,.k-5)+ Ay (is +sk-5) +B (ii+e4+5)4+ Fh (L7-5,k +5) (22) 
2HB (i,j,k) = Ay(ij+s,k) +A (ij-3,k) (23) 


- Ey(i,j+sk+1)-Ey(ijtsk)  Ey(ij+—k+1)-E,(ij,k+— 

HBG j,k) = A, (ii +5,k) +8 (eee —2 zs Aiea) al nay a ) u 2) (24) 
2.3. Shielding Effectiveness by MOM Method: 

The curved section in this enclosure needs accurate analysis and simple numerical methods will not 
be responsive. The MOM numerical solution method can solve this problem. The way of solving in this 
method is the basic equations of electromagnetic fields, which fields are obtained in a boundary condition by 
numerical solution of these equations. Finally, solving the matrix (depending on the type of problem and its 
convergence) is done by numerical solution. To solve the divergence problem of the integral expansion at 
high frequencies, the combination of physical optics and MOM is used in simulation. This simulation has 
been performed at 0-3 gigahertz of frequency range so we can say that the combination of MOM and 
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physical optics will respond to the simulation of electromagnetic fields in complex structure at high 
frequencies [17]. The field equation can be written as Equation 25 whit considering the existence of 
Ei(—M,) + Eiji) =0 in the enclosure and EQ(ji™® ) + E°(M,) + E°(j°) =0 out of the enclosure. 
Helmholtz's equation for correction of electric field is used in two parts of rotational with a zero divergence 
(E,,) and non- rotational with a zero Curl (F,,). The magnetic field and total electric field are calculated and 
represented with using the Helmholtz equation. The non-circular part is empty, so the study of the circulation 
section is sufficient. The basic idea of MOM is that the unknown quantity f is expanded to a set of Linear 
Independent Functions which is approximated by the following limited series (Equation 26). a, 1s the weight 
factor. After considering the boundary conditions, we define the internal multiplication or a moment between 
the fundamental function f, (1r’) and the weight function f,,(r) which is shown in Equation 27 [12, 16]. 


He U'™’ ) + Hp (Ms) + He U°) = Hi(—-Ms) + HU") (25) 
f oe vei Antn (26) 
(fm fn) = Sp fm) Sp fn dr’ dr (27) 


3. SIMULATION AND COMPUTING 

Various magnetic materials such as iron, copper and zinc are commonly used for equipment as 
passive shielding. The accurate use of the shielding method can be used as a way to reduce and eliminate the 
destructive effects of electromagnetic fields on equipment. Shielding can be placed around circuits, systems, 
wires or communication cables [23]. The reduction of the radiant electromagnetic field in the shielding is 
determined by the SE coefficient [24-27]. An enclosure with dimensions of 20*20*40 is simulated with the 
aperture and the simulation of Shielding effectiveness is performed by using three methods of numerical 
solution, FDTD, MOM and equivalent circuit in a Curved enclosure with Apertures. This simulation is 
performed in four steps in the enclosure. One time, this simulation was performed with a large 10*20 cm 
wide aperture (Figure 5) and then with 12, 36 and 60 small apertures (Figure 7-9). The SE curve for the 
magnetic and electric field is shown in Figure 6. Finally, the values of different numerical solving methods 
are compared with the measured values and are shown in Figure 10. 
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Figure 6. SE curve for (a) the magnetic field (b) the electric field 
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Figure 7. Enclosure with 12 apertures (2*2 cm) (b) SE curve 





N 


SE-Efieid 
— MOM Method 
--| ~ FDTD method 
— Circuit Model 








SE-Efeld(dB) 











re ow 100 
SE-Eftela 
— MOM Method 
| ~ FDTD method 
~ Circuit Model 


o 
SE-Efield( dB) 





00 05 10 15 20 25 3.0 








100 = 
SE-Efield 
— MOM Method 
~ FDTD method 
~ Circuit Model 
| — Measured 


SE-Efieid(dB) 








40 
0.0 05 1.0 15 20 25 3.0 


Figure 10. Comparison between three numerical solution methods with measured values of SE curve with 60 
apertures (2*2 cm) 
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In the final section, we attempt to evaluate the SE with using different materials and measure their 
impact on the quality of the enclosure efficiency. We have chosen three materials, Iron, Copper and 
Mu-Metal. Simulation of SE curves with different materials is shown in Figures 11-13. 
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Figure 11. SE E and H field curves in the Iron enclosure 
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Figure 12. SE E and H field curves in the Mu-Metal enclosure 
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Figure 13. SE E and H field curves in the Copper enclosure 


4. CONCLUSION 

In this paper, the effect of different apertures and material on the shielding effectiveness in a 
complex structure has been accomplished by three methods of numerical solution, FDTD, MOM and 
equivalent circuit. The complex structure of the enclosure, despite its curvature in the enclosure, allows the 
equipment to be minimized and optimizes the dimensions. in fact, The use of apertures in the wall of the 
enclosure is due to the modeling of the connection of cables and ventilation. The enclosure simulation was 
performed to determine the shielding effectiveness in four different compartments and three different 
materials in frequency range of 0-3 gigahertz. What is clearly visible is the high level of effectiveness at low 
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frequencies in all three numerical methods. Increasing the number of apertures and decreasing its dimensions, 
makes the shielding better, but the number of resonance points also increases. The simulation shows that the 
MOM method is closer to the measured values for most of the selected frequency range, which explains on 
the accuracy of this method in complex structures. 

The use of different materials in the production of the enclosure shows that the Mu-metal enclosure 
is better and the SE in this enclosure is optimally increased throughout the entire frequency range. The 
presence of resonant points in frequency range, The importance of choosing the exact frequency of 
equipment is shown. 
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